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ABSTRACT 

Ground-based long-term optic variability of AM Her, covering the period 
between 2003-2008, has been conducted to study the features seen in both low 
and high states of the system. Low-state analysis shows the presence of short- 
term, low-amplitude light variations of about 0.02-0.03 mag with a mean power 
time between 16 s and 226 s. Brightness variations on the order of 0.7-2 mag, 
which could be due to the stellar activity of the component in the system, are 
also detected. A total of 30 years times of minimum light given in the literature 
are combined with nine times of minima obtained in this study. We represented 
the (observed-calculated) diagram by a parabolic curve and also by two broken 
lines. Under the assumption of a parabolic variation, we estimate an increase in 
period, ^ = 7.5(1.2) x 10^^ days yr~-^, with a mass transfer rate of M = 8(2) x 
lO~^M0?/r~l, in agreement with the previous findings by a different method. 

Subject headings: Binaries: close — binaries — stars: Cataclysmic Variables — 
stars: individual (AM Her)- techniques: photometric 



^Bascd on observations gathered with the 1.5m Russian-Turkish telescope (RTT150) and ROTSE Illd 
at the TUBiTAK National Observatory. 
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1. Introduction 

AM Her systems (polars) are semidetached binaries that consist of strongly magnetic 
white dwarf (WD) primaries and red dwarf (RD) secondaries. Polars were first recognized 
in 1976 with the discovery of circular polarization in AM Her (Tapia 1976, 1977). Magnetic 
fields play a crucial role in determining the system's parameters. The field strength of the 
primary is so high that the material fiowing from the companion does not form an accretion 
disk around the WD, but is guided along the field fines to an accretion column that forms 
near the magnetic pole of the primary. The flux distribution from the column consists of 
hard X-ray bremsstrahlung, an approximately blackbody spectrum in the UV and soft X-ray 
and cyclotron emission, which is the primary source of optical radiation. 

The brightest polar AM Hercufis (RX J1816.2 +4952 = EUVE J1816+49.8 = 3U 
1809-1-50 = H 1816-1-49) was classified as a cataclysmic variable (CV) by Berg & Duthie 
(1977). It is not an eclipsing binary despite the observed large- amplitude minima in the 
light curves. The system has been studied for different wavelengths. Each study revealed 
different characteristics of the polar. The system shows long-term, non-periodic variations 
where the brightness of the system varies by about 2 mag (see, Hessman, 2000), known as 
high and low states. High and low states of polars are thought to be due to the variation in 
the mass transfer rate from the RD to the WD. In a recent study on the high and low states 
of the system by Wu & Kiss (2008) it was found that the magnetic field of the primary is a 
crucial parameter in regulating these states. Previously, Livio & Pringle (1994) explained the 
observed low state of the system with the starspots migrating under the inner Lagrangian 
point (LI). 

Observational properties of polars generally depend on the observed state and wave- 
length. For instance, short-term, low-amplitude variations in the X-ray/optical bands, known 
as ffickering, are detected when the system is in high and intermediate states (King, 1989; 
Bonnet-Bidaud et al. 1991). Different mechanisms have been suggested to explain these 
variations. Szkody & Margon (1980), using the cross-correlation functions of high state ob- 
servations of AM Her, reported the strong correlation in Johnson U, V, and A4686 features 
and discussed the ionizing radiation as a responsible mechanism for the observed small- 
amplitude brightness variations. Another mechanism is the oscillation of the magnetic fiux 
tubes (Tuohy et al. 1981). Larsson (1988), on the other hand, proposed an oscillatory shock 
height model to explain optical variations with periods, of a few seconds. King (1989) dis- 
cussed that X-ray irradiation of the accretion flow below the LI point produces oscillating 
ionization fronts. These ionization fronts modulate the accretion rate through LI. The 
timescale of these oscillations is the dynamical timescale near the LI, point which is about 
8 minutes for AM Her. Besides the high state, the low state of the system has also been the 
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subject of interest because of its physical properties and the poorly defined characteristics 
that can change over time. 

Complex and unpredictable observational properties of polars prompted us to obtain 
long-term optical variation of AM Her. Studying the light variations obtained over a long 
period of time is necessary to see the complex structure. The study of AM Her was carried 
out with the 1.5m RTT150 and ROTSE Hid (Robotic Optical Transient Search Experiment- 
Hid) telescopes of the TUBITAK National Observatory (TUG). The results are presented 
in Section 2. In Section 3 we analyze the observations and discuss them in Section 4. 



2. New observations and light variation 

Optical photometry of the system was obtained using the Russian-Turkish 1.5m tele- 
scope (RTT150) over 19 nights between the period 2003-2007 (Tabled]). All the images were 
obtained using the Andor CCD. The Andor CCD camera is equipped with a set of Cousins 
(Rc, Ic) and Johnson {V) filters. During the data reduction a few comparison and check stars 
were chosen in the same CCD frame including CSC 3533 1026 and CSC 3533 1021. All CCD 
reductions were done with the IRAI^ package. In Fig. [1] the times of RTT150 and ROTSE 
Illd observations are indicated by vertical lines with the optical light curve of AM Her from 
the data of the American Association of Var,able Star Observers (AAVSO). Light variations 
of AM Her, obtained with RTT150, over long periods of time show that the light- variation 
amplitude changes over time. The light curves, in Fig. [21 are plotted as a function of the 
Julian Date (JD) to make clear any possible variation between successive orbital phases (see, 
Kalomeni et al. 2005 for the light variation of the system obtained in 2003). 

ROTSE Ilici was used to carry out long-term variation of the system (Tabled]). Fig. [3^ 
shows long-term light variation of AM Her spread over 205 nights, obtained with the ROTSE 
Illd, between the period 2005-2008. During the low state two noticeable brightness varia- 
tions have been detected. Low-state observations of AM Her exhibit rather weak brightness 
variations with respect to the high and intermediate states. During the low-state the mass 
transfer from the secondary is thought to decrease or cease. If the accretion is almost negligi- 
ble then the characteristic features of the component can exhibit themselves in observed light 
curves (Kafka et al. 2005). Therefore, because the secondary is a late-type main-sequence 



^IRAF is distributed by National Optical Astronomy Observatories, which is operated by the Association 
of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation, U.S.A. 

^The ROTSE III system's details described in Akerlof et al. (2003). 
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Fig. 1. — Long (more than 25 years) optical light curve of AM Her from AAVSO data. The 
times of RTT150 and ROTSE Illd are shown with vertical lines (see Table 1 for the data) 
. (A color version of this figure is available in the online journal.) 



Table 1: Summary of the observations of AM Her with RTT150 and ROTSE Hid* telescopes. 
HJD+ shows JD start+2400000 



Run 


HJD+ 


Filter 


Nobs 


State 


1 


52858.46-52858.59 


Rc 


162 


Intermediate 


2 


52859.47-52859.59 


Rc 


195 


Intermediate 


3 


53037.55-53037.67 


Rc 


148 


Active 


4 


53193.29-53193.57 


Rc 


444 


Low 


5 


53194.30-53194.52 


Rc 


261 


Active 


6 


53195.33-53195.59 


Rc 


371 


Low 


7 


53425.51-53425.55 


Rc 


61 


Low 


8 


53426.47-53426.53 


Rc 


70 


Low 


9 


53450.53-53450.63 


Rc 


129 


Low 


10* 


53464.95-53582.84 




387 


Low + active + high 


11 


53489.45-53489.60 


Rc 


229 


Active 


12 


53682.20-53682.33 


Rc 


90 


Low 


13 


53682.20-53682.30 


V 


51 


Low 


14 


53683.19-53683.35 


Rc 


266 


Active 


15 


54014.23-54014.31 


Rc 


138 


Active 


16 


54015.26-54015.40 


Rc 


200 


Low 


17 


54102.66-54102.68 


Rc 


14 


Low 


18 


54103.17-54103.20 


Rc 


16 


Low 


19 


54324.40-54324.56 


Ic 


253 


High 


20 


54324.40-54324.60 


Rc 


297 


High 


21* 


54012.22-54402.21 




1020 


Low + active -|- high 


22* 


54089.21-54573.53 




70 


High 
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Light curves of AM Her obtained between the period 2004-2007. All hght curves 
except (O) that is obtained in Ic. 
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star, we can expect to detect stellar activity in the low state. Such a brightness variation in 
AM Her was observed by Shakhovskoy et al. (1993). They reported an approximately 2 mag 
flare event with a 20 minutes duration (see also, Bonnet-Bidaud et al. 2000). We detected 
similar events of 2.14 mag for 1 hour and 20 minutes in 2005 with ROTSE Illd (Fig. [Sb)- 
Following this variation ROTSE Illd detected a 1.37 mag excess in about 30 minutes on 
July 20, 2005 (see. Table [2]). During AM Her's low-state relatively small amplitude flaring 
with amplitudes of 0.2-0.6 mag and lasting 15-90 min was reported by Kafka et al. (2005). 
Likewise, during the observing runs performed between 2004 and 2006 similar variations 
were detected (Table |2l Fig. H]). Magnitude excesses, owing to the possible flare events, with 
respect to the quiescence level are shown in Table [2j In the RTT150 observations we have 
also determined times of minima, derived by the Kwee-van Woerden method (Kwee & van 
Woerden, 1956) and for the asymmetric minima by freehand curve, from the individual light 
curves. In Table E] the times of minima are shown. The errors in times of minima are of the 
order of 0^.0002 - 0^.0007. 



3. Period Change Analysis 

Times of minima in AM Her, as well as in other polars, is known to offset (e.g. Bailey 
et al. 1993). The nature of the observed shift in times is poorly understood. However, if the 
orbital period of the system is determined accurately and the WD is synchronized with the 
orbital period, then these shifts generally attributed to the oscillation of the magnetic pole 
(Bailey & Axon, 1981; Bailey et al. 1993). Any variation in the mass accretion rate alters 
the accretion geometry. In this case while the location of the pole remains flxed the position 
of the spot with respect to the pole changes (Cropper 1989). On the other hand, the orbital 
periods of interacting binaries are known to change because of different processes. One of 
these is the mass transfer between the components. In AM Her systems, the RD component 
loses mass to the primary WD star. In semidetached binary systems, the displacement in 



Table 2: Magnitude variations for AM Herculis. HJD* shows JD start +2400000 



HJD* 


Filter 


Phase 


magnitude excess 


53037.5712 


Rc 


0.1 


0.7 


53037.6200 


Rc 


0.5 


0.71 


53194.4846 


Rc 


0.3 


0.12 


53489.5381 


Rc 


0.8 


0.5 


533545.8400 




0.52 


2.14 


53571.8293 




0.01 


1.37 


53683.2006 


Rc 


0.8 


0.45 


53683.2484 


Rc 


0.2 


0.7 


54014.2494 


Rc 


0.56 


0.7 
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minima times causes a parabolic variation, either upward or downward, in the difference 
between the observed (O) and calculated (C) time of minima diagram. Thus, analysis of the 
observed times of minima is important to determine any variation, due to mass transfer, in 
the orbital period. Unfortunately, there are almost no studies done on (0-C) variation of 
other polars in the literature. Previous (0-C) studies of AM Her were performed by Young 
& Schneider (1979) and Mazeh et al. (1986). The study of Young & Schneider (1979) shows 
no evidence for any continuous period variation. The latter study by Mazeh et al. (1986) 
shows a downward curved parabola with P/P=—5 x lO^^'^s^^ We collated the additional 
minima times obtained since then with those obtained in this study (Table [3]) to revise the 
(0-C) variation. We assigned the same weight for all minima points during analysis. The 
starting epoch for the primary minimum was adopted from the Szkody & Brownlee, (1977). 
The (0-C) diagram of AM Her constructed with an initial light element can be represented 
by the relation, 

HJD Mini = 24 443014.7136(2) + 0.128927048(2) x E 

+ 1.33(20) X 10-12 xE^. ^ ' 

The (0-C) variation can be represented both by a parabola and with two broken lines. 
First, if the long-term variation fits to a parabola, as is expected for binaries where mass 
transfer take place, the solid curve in Fig. [5^ is for the best solution to the (0-C) variations. 
Then, the resulting solution with a quadratic term is an indication of the effect of mass 
transfer in AM Her. The mass transfer in polars is from the less massive secondary to the 
more massive primary WD star. This increases the orbital period in the conservative case 
and the resulting (0-C) curve is a parabola with a positive quadratic term. On the other 
hand, as a result of the nonconservative mass transfer and magnetic activity of the late- 
type component, mass loss from the system may also occur in polars. However, the upward 
curving parabolic variation in the (0-C) diagram of AM Her indicates that the dominant 
effect is the mass transfer from the RD to the WD. If this variation fits to a parabola then the 
corresponding rate of period increase is dP/dt = 7.5(1.2) x 10"^ days yr-^ with a conservable 
mass transfer rate of M = P/3P[MiM2/(Mi - Ms)] = 7.6(2.3) x lO'^M© jt-\ for a WD 
with a mass O.88M0 (Bailey et al. 1988), which agrees with the maximum mass transfer rate 
given by Hessman et al. (2000) assuming that this variation is caused by stellar spots. On 
the other hand, the P/ P value is of the order of another recently studied polar (Andronov & 
Baklanov, 2007). We can also fit the data with two broken lines. One of them is horizontal 
while the other one indicates an increase in the period. Nevertheless, observations in the 
next decades are necessary to clarify the true shape of the (0-C) diagram. 

The accretion luminosity of the WD is L = —GMiM/Ri where Mi is the WD mass, 
i?i is the radius of the WD, and M is the mass accretion rate. Using the M estimated from 
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the possible parabolic (0-C) variation, we can determine the accretion luminosity to derive 
the Alfven radius. The Alfven radius is given by (Prank et al. 2002) 

T, = 2.9 X lQ^Ml/'Rf"L-,^"B%\ (2) 

where is the radius of the WD in units of lO^cm, L37 its luminosity in units of 10^^ erg 
s~^ and B12 is the surface magnetic field strength in unit of 10^^ G. The estimated Alfven 
radius for the adopted value of M = O.88M0 is 2.02 x 10^°cm . 

AM Her shows evidence for brightness variations on a timescale of minutes (approx- 
imately 4.5 minutes; e.g., Bonnet-Bidaud et al. 1991). One of the explanations for the 
origin of the quasiperiodic X-ray variations observed in AM Her is discussed by Tuohy et 
al. (1981). They presented the oscillation of the magnetic flux tubes, through which matter 
flows to the stellar surface, as a possible mechanism for these variations. If these oscillations 
occur at the point where the matter channeled by the magnetic field, they can lead to the 
quasi-periodic variations in the accretion rate. The timescale for an Alfven wave to cross 
the magnetosphere is 

P„,e(r) = 2 X lQ-V^''Ly:rTM;"'RfB^' s, (3) 

where rg is the radius where the flowing matter is channeled and quasiperiodic variations 
may arise, / is the fraction of stellar surface where the accreting matter flows. Adopting 
/ = 9.1 X 10~^ as the characteristic dimensionless size of the flow (Hessman, 2000), we flnd 
the related r on the timescales of interest as 2.03 x 10^°cm. The r estimated for the 4 minutes 
oscillations agrees very well with the estimated Alfven radius. This indicates that this model 
adequately describes the observed brightness variations. 

The Lagrangian radius of the WD is Rli/cl — 0.5 — 0.227 log^ (Plavec & Kratochvil, 
1964) with q the mass ratio of the components and a the orbital separation. For q — 0.31 
and a = 7.85 x 10^°cm, this yields Rli = 4.8 x 10^°cm and f» 0A2Rli. Therefore, as is 
expected both the magnetospheric radius and the radius where the matter channeled toward 
the WD are smaller than the Lagrangian radius (Ferrario et al. 1989). The results obtained 
in this study are in agreement with the results presented by Bonnet-Bidaud et al. (1991). 
They discussed that their observed 270 s variations correspond to a radius ~ 2.1 x 10^°cm 
for AM Her. The flickering timescale depends on M and / since Mi, R and B cannot 
change on the timescales of 3-8 minutes. Hence, as Bonnet-Bidaud et al. (1991) reported, 
any inhomogeneities in the accretion matter can produce these brightness variations on the 
timescales of interest. 



- 9 - 



4. Results and Discussion 

AM Her-type systems show large-amplitude variations over years (Fig. 1-3) as well as 
short-term low-amplitude variations. In this study, five years of observations obtained in both 
states of the system are presented. Low-state photometry reveals weak orbital modulation 
but occasionally fiaring-type variability of the secondary. Large fiare events, as well as 
smaller amplitude fiares, are detected. Three of the nine magnitude excess events detected 
fell within the primary minimum and one within the secondary minimum. Durations of 
brightness variations range from tens of minutes to an hour. Flickering with an amplitude 
of 0.01-0.60 mag occurs on a timescale of at least a few minutes. 

We have obtained a total of nine times of minima, using them we could perform a 
period analysis of the system. The times of minima in AM Her show shift, which is, generally, 
assumed due to obscuration of the post-shock radiation of the main accretion column (Bailey 
et al. 1993). These changes are thought to be responsible for the observed shift in (0-C) 
diagram. On the other hand, AM Her systems are described with a mass-donating RD and 
a magnetic WD star, therefore, they are classified as semidetached binaries. Semi-detached 
binaries known to show a parabolic variation in the (0-C) diagram due to the mass transfer 
between the components. Using the available times of minima, we find for the first time, an 
evidence for an upward parabola general property observed in semi-detached binaries where 
the matter fiows from less massive to the more massive. Using this variation we derive an 
orbital period evolution time of about 1.7x10^ yr. In addition, the (0-C) variation can be 
described with two broken lines. A total of 30 years times of minima are collated. This time 
is quite long to see any variation in the (0-C) diagram. However, it is apparent from Fig. [5^ 
that the upward parabola is not so clear as seen in binaries with nondegenerate components 
(e.g., Kalomeni et al. 2007). If the period variation is due to the conservative mass transfer, 
the mass transfer rate between the components is M = 7.6 x 10~^ MQyr~^. This mass 
transfer rate is small in comparison with the binaries with non-degenerate components (ibid). 
Using timescales for the gravitational radiation and magnetic braking (Yakut et al., 2008) 
we have calculated the gravitational radiation timescale and magnetic braking timescale of 
AM Her about 7 Gyr and 1.3 Gyr, respectively. The mass accretion rate timescale of the 
system is much less than these. Thus the gravitational radiation is not as important as the 
other mechanisms for AM Her. On the other hand, if the orbital period of AM Her were 
half of its present value then gravitational radiation and magnetic braking would be much 
more important. At periods of Ih, 2h, and 8h (e.g., V1309 Ori) the gravitational radiation 
timescales are about 0.4 Gyr, 2.4 Gyr, and 96 Gyr, respectively. In AM Hers, the binary 
geometry, period change ratio, angular momentum loss mechanisms, etc. can all change 
because of the strong magnetic field of the WD and the interaction between components' 
magnetic field (see also Wickramasinghe & Wu, 1994; Webbink & Wickramasinghe, 2002). 



- 10 - 




Fig. 3. — (a) The figure shows long term hght variation of AM Her spread over two hundred 
five nights obtained with the ROTSE Illd between the period 2005-2008, (b) the magnitude 
excess observed in HJD2453545 



Table 3: The times of minima of AM Her in HJD* (HJD - 2400 000). 



HJD* Passband Ref. HJD* Passband Ref. 



43014.71266 


V 


1 


44133.02540 


V 


4 


43014.84127 


V 


1 


45591.32259 


V 


5 


43015.74554 


V 


1 


45600.34559 


V 


5 


43015.87731 


V 


1 


46000.2788 


V 


5 


43031.72862 


V 


1 


46001.31151 


V 


5 


43031.86055 


V 


1 


46132.55800 


V 


5 


43032.6336 


I 


2 


51277.5182 




6 


43033.661 


I 


2 


51708.45993 


980-1180A 


7 


43062.5439 


I 


2 


52858.5548 


Rc 


8 


43062.8024 


V 


1 


52859.5259 


Rc 


8 


43069.6354 


I 


2 


53193.3160 


Rc 


8 


43083.5591 


I 


2 


53193.5168 


Rc 


8 


43635.88762 


6900-7400A 


3 


53195.4500 


Rc 


8 


43636.78853 


8200-8700A 


3 


53450.5989 


Rc 


8 


43704.73386 


8000-8350A 


3 


53489.5276 


Rc 


8 


43704.86349 


8000-8350A 


3 


53683.2436 


Rc 


8 


43705.89342 


7500-7800A 


3 


54015.2955 


Rc 


8 



References for Table [31 1-Szkody & Brownlee (1977), 2-Olson (1977), 3- Young & Schneider 
(1979), 4- Young et al. (1981) based on Bailey & Axon (1981) observations, 5-Mazeh et al. 
(1986), 6-Safar & Zejda (2002), 7-Hutchings et al. 2002 (the average minimum time is used), 
8-This study. 
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Fig. 4. — Expanded views of the Rc- band brightness variation during the low-state of AM 
Her (a) April 28, 2005, (b)-(c) November 8, 2005, and (d) October 5, 2006. 
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Fig. 5. — a) Observed-calculated times of minima of AM Her vs. epoch, and b) the residuals 
of parabolic variation vs. epoch (see text for details). (A color version of this figure is 
available in the online journal.) 
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